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PRE-SEQUENCES OF 
MATRIX ORTHOGONAL POLYNOMIALS 

JUAN TIRAO 


Abstract. We introduce the notion of a pre-sequence of matrix orthogonal 
polynomials to mean a sequence { A, } r ,>o of matrix orthogonal functions with 
respect to a weight function W, satisfying a three term recursion relation and 
such that det(Fo) is not zero almost everywhere. By now there is a uniform 
construction of such sequences from irreducible spherical functions of some 
fixed A'-types associated to compact symmetric pairs (G, K ) of rank one. Our 
main result is that { Q n = A„Fq -1 }„>o is a sequence of matrix orthogonal 
polynomials with respect to the weight function FqWFq, see Theorem 12.II 


1. INTRODUCTION 

In the seminal paper [5] a sequence of matrix polynomials {F n } n >o was con¬ 
structed from the irreducible spherical functions of the pair (G', K) = (SU(3), U(2)), 
for each fixed X-type. More precisely, the restrictions of all these functions to the 
domain of an appropriate local coordinate x of an Iwasawa subgroup A of G, are 
carefully arranged as the columns of the matrix polynomials F n (x) , 0 < x < 1, 
n > 0. The irreducible spherical functions of a pair (G, K), where G is a compact 
group and I\ a closed subgroup, of a fixed X-type tt are in a one to one corre¬ 
spondence with the unitary irreducible representations of G that contain tt as a 
subrepresentation. In our case, from Schur orthogonality relations of unitary ir¬ 
reducible representations of SU(3) it follows that is a set of orthogonal 

polynomials with respect to a weight function W = W(x) supported on the closed 
interval [0,1]. Besides, the irreducible spherical functions are eigenfunctions of the 
algebra D(G) K of all left invariant differential operators on G which are also right 
invariant under K. In particular, this implies that ( F n D){x ) = A n F n (x), for all 
0 < x < 1, where D is a second order right-hand side linear differential operator 
with matrix valued coefficients. Moreover, using the explicit branching rules for 
the tensor product of the fundamental spherical representation and an irreducible 
unitary representation of G, one proves that each sequence {F n } n > o satisfies a three 
term recurrence relation of the form xF n = A n F n -\ + B n F n + C n F n+ 1 , where the 
coefficient matrices are explicitly known. Such sequences {X„}„>o are very close to 
being sequences of matrix orthogonal polynomials, but they are not, because F n is 
not of degree n and its leading coefficient fails to be nonsingular. 

The results in mum yield for the first time examples of matrix valued classical 
orthogonal polynomials {$l a ’^}„>0 of Jacobi type, where the parameters a, /3 take 
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the values a £ Z>o and /3 = 1; classical in the sense that they are eigenfunctions 
of a fixed second order differential operator. This was attained from the sequences 
{Fn}n> o> mentioned above, by defining the matrix valued function by means 
of = F u Fq 1 . This occurred after a useful conversation of F. A. Griinbaum 
with A. Duran that steered him to this definition, see Section 5 of [7]. These were 
the first examples of such classical sequences of matrix orthogonal polynomials 
which were sought for by the experts since the publication [lj of A. Duran. Then 
in [4] F. A. Griinbaum gave a completely explicit description of an extension of 
these sequences to arbitrary values of the parameters a,/3 > —1, in the 

case of size two. The sequences {'F[, Q ' d) } n >o become classical sequences of matrix 
orthogonal polynomials. After that, many more examples of these of polynomials 
were found by analyzing different compact symmetric pairs of rank one, see for 
example ng ezi uni na mi na ng tini nsi na], and by some ad-hoc methods as in 

m. 

In 120] through a trial and error process we found a function ^ used to play the 
same role as Fq in i 13 El i, in describing the irreducible spherical functions of 
any AT-type associated to the symmetric pairs (SU(3),U(2)) and (SU(2,1),U(2)). 
The same function T was used in [□HE to build sequences of classical orthogonal 
polynomials from sequences of matrix polynomials {F n } n >q constructed from irre¬ 
ducible spherical functions of the pairs (SU (d + l),U(d)). This same strategy was 
used in [iacsici]. 

The sequences of matrix polynomials {Fnln^o constructed from the irreducible 
spherical functions of compact symmetric pairs (G, K ) of rank one, for some fixed 
A'-types, are sequences of matrix orthogonal functions with respect to a weight 
function W, satisfying a three term recursion relation and such that det(F 0 ) is not 
zero almost everywhere. Under these conditions we prove that {Q n = A , T jF 0 _1 } n >o 
is a sequence of matrix orthogonal polynomials with respect to the weight function 
F 0 WFq, see Theorem O For this reason such sequences are called pre-sequences 
of matrix orthogonal polynomials. 

The simple and fundamental fact expressed in Theorem 12.11 which is the core 
of this paper, was conceived in January 2012 and immediately shared with my 
coauthors. At the time Maarten van Pruijssen was working on his Ph.D. thesis 
where he extends the results obtained by many of us for (SU (d + l),U(d)) and 
(SO (d + l),SO(d)), to other compact Gelfand pairs (G, AT) of rank one. But in 
these new cases he needs to restrict his attention to those A'-types 7r such that 
ind^-(7r) is a multiplicity free direct sum of irreducible representations of G. In this 
very nice extension he uses our idea in defining the function T = Fq, see Theorem 
3.5.12 in 13(, also consult [10, 14] , 


2. Pre-sequences of MOP 

By a weight matrix of size N on the real line W = W (x) we mean an integrable 
function on an interval (a, b ), such that W ( x ) is a (self-adjoint) positive semidefinite 
N x Af-matrix for all x £ ( a,b ), which is positive definite almost everywhere and 
with finite moments of all orders: for all n £ No we have 


x n W{x)dx £ MatAr(C). 
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More generally we could assume that W is a Borel measure on the real line of 
N x IV-matrices, such that: W(X) G Mjv(C) is positive semidefmite for any Borel 
set X, W has finite moments of any order, and W is nondegenerate, that is for P 
in the polynomial ring Matjv(C)[x] 

(P,P)= [ P(x) dW(x)P(x)* = 0, 

it 

only when P = 0. If M G Mat/v(C), M* denotes the conjugate transpose of M . 

Given a weight matrix W = W ( x ) of size N on the real line, supported on 
the closed interval [a, &], we introduce as in [IT] and jl2| the following Hermitian 
sesquilinear form in the linear space Mat jv (C) [x], 

(P,Q)= f P(x)W(x)Q(x)* dx. 

J a 

The following properties hold: for all P,Q,R G Matjv[x], T G Matjv(C) and a, b G C 
we have 


(aP + bQ , R) = a(P , R) + b{Q , R), 

(TP,Q) =T(P,Q), 

(P,Q)*=(Q,P), 

(P, P) > 0; if (P, P) = 0 then P = 0. 

These properties follow directly from the definition except the last one (cf. Proposi¬ 
tion 2.2 in M)- In other words Matjv(C)[x] is a left inner product Matjv(C)-module. 

By a sec^uence of matrix orthogonal polynomials {Q n }n> o of size N we mean a 
sequence of matrix polynomials such that deg(Qn) = n, the leading coefficient of 
Q n is nonsingular and ( Q m , Qn) = 0 for all m ^ n. Then one establishes that there 
exists a unique sequence of monic orthogonal polynomials {P n }n>o and that any 
sequence {Q n }n> o of matrix orthogonal polynomials is of the form Q n = M n P n 
where M n G GLjv(C) is arbitrary for each n > 0 (cf. Corollary 2.5 in j9j). 

Two matrix weights are said to be equivalent W' ~ W, if W' = MWM* for some 
M G GLat(C). In such a case the map {Q n }n> o ^ o establishes a 

bijection between the set of all sequences of matrix orthogonal polynomials with 
respect to W and the corresponding set with respect to W'. 

A standard argument, given for instance in mm shows that a sequence of 
matrix orthogonal polynomials {Q n }n> o satisfies a three term recursion relation 

(1) xQ n (x) — A n t^ n _i(x) -f- T C n Q n -\. i(x), n ^ 0 

where Aq = 0 and C n is a nonsingular matrix . 


We will consider a sequence {F n }n>o, of matrix valued functions on (a, b) such 
that det(Po)(x) is not zero almost everywhere. Let W = W(x) be a weight matrix 
on (a, b) except that we not require W to have finite moments of all orders. If 
besides Fi(x)W(x)Fj(x)* is an integrable function on [a, b\ for all i, j > 0 and 


Fi(x)W(x)Fj (x)* dx = 0 for all 


*7 1 3, 


such a sequence will be called a sequence of matrix orthogonal functions. 
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Theorem 2.1. Let {F n } n >o be a sequence of matrix orthogonal functions on ( a,b ) 
with respect to a weight function W supported in [a, b}. If the sequence satisfies a 
three term recursion relation of the form 

xF n = A n F n _i + B n F n + C n F n+ i n > 0 

with Aq = 0 and C n nonsingular, then F n = Q n Fo where Q n is a uniquely deter¬ 
mined polynomial of degree n with a nonsingular leading term, for all n > 0. The 
moments ofW' = FqWFq are all finite and the sequence {Q n }n >o is a sequence of 
matrix orthogonal polynomials with respect to W', satisfying 

(2) xQn = AjiQji —i -f- B n Qn -(- CnQn -\-1 n ^ 0. 

Proof. By induction on n > 0 we will prove that F n = Q u Fq where Q n is a 
polynomial. For n = 0 we have Fq = QqFq which is true for Qq = I. Now we 
assume that F v = Q v Fq with Q v polynomial for all 0 < v < n. Then we have 

xQ n Fo — xF n = A n F n —i B n F n + C n F n +i 
= A n Q n _iF[j ~h B n Q n F 0 T C n F n _ j_i. 

Now by hypothesis C n is nonsingular for all n > 0. Then we have 

l'h ■ 1 — c n {xQ n A n Cf n —i BnQrf) Fq. 

Thus we can take 

(d) Qn+1 — C n (xQ n A n (f n — x B n Qrfj, 

completing the proof of the inductive step. Now from ([3]) it follows by induction on 
n > 0 that Q n is a polynomial of degree n with a nonsingular leading coefficient. 

If a matrix polynomial P = P{x) satisfies P{x)Fq{x) = 0 for all x, then P(x) = 0 
for almost all x, because by hypothesis Fq(x) is not singular for almost all x £ (a, b). 
Hence P = 0. Taking this into account, the three term recursion relation satisfied 
by {Qn} n >o is obtained from the one in the hypothesis, by canceling Fq from both 
side of such equations. Similarly the uniqueness of the Q n follows from the same 
reason. 

Moreover, for all i,j > 0 we have 

/ Qi(x)W'(x)Qj(x)* dx = / Qi(x)F 0 (x)W(x)F 0 (x)*Qj(x)* dx 
J a J a 

= f Fi{x)W{x)Fj(x)* dx = 5ij{F i ,F j ). 

J a 

Hence {Q n }n> o is a sequence of matrix orthogonal polynomials with respect to W' 
which implies that all its moments are finite. The theorem is proved. □ 

It is worth to observe that the hypothesis of integrability of Fi(x)W{x)Fj(x)* 
for all i,j > 0 can be replaced by the hypothesis that F 0 (x)W(x)Fq(x)* has all its 
moments finite. 

We feel appropriate to refer to a sequence {Fkj^o of matrix orthogonal functions 
that satisfies a three term recursion relation as m as a pre-sequence of matrix 
orthogonal polynomials. 

We come now to the notion of a differential operator with matrix coefficients 
acting on matrix valued polynomials, i.e. elements of Matjv(C)[a;]. These operators 
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could be made to act on our functions either on the left or on the right. One finds 
a discussion of these two actions in [TJ. The conclusion there is that if one wants 
to have matrix weights that do not reduce to scalar weights and that have matrix 
polynomials as their eigenfunctions, one should settle for right-hand-side differential 
operators. We agree now to say that D given by 


D = Z20 i F i ( x ), 

i =0 



acts on Q(x) by means of 


QD = Y J d\Q){x)F i (x). 

*=o 

Given a sequence of matrix orthogonal polynomials {Q n }n> o with respect to a 
weight matrix W = W{x), we introduced in [9] the algebra D(W) of all right-hand 
side differential operators with matrix valued coefficients that have the polynomials 
Q w as their eigenfunctions. Thus 


(4) V(W) = {D : Q n D = T n {D)Q n , T„(D) € Matjv(C) for all n > 0}. 


The definition of T>(W) depends only on the weight matrix W and not on the 
sequence {Q n }n> o, because Q n = M n P n for some M n £ GLat(C) where {P n } n > o 
is the sequence of monic orthogonal polynomials, besides , if D £ T>(W), then D 
is hypergeometric, meaning that the coefficient of the jth-order term is a matrix 
polynomial of degree less or equal to j, see Proposition 2.6 of [9j. 

If W' ~ W, say W' = MWM *, M £ GL Ar (C), then the map D ^ MDM~ X 
establishes an isomorphism between the algebras V(W) and V(W'). In fact, if 
{Qn} n >Q is a sequence of matrix orthogonal polynomials with respect to W, then 
{Q' n = MQ n M~ 1 } n > o is a sequence of matrix orthogonal polynomials with respect 
to W'. Moreover, if Q n D = T n (D)Q n , then 

Q'JMDM- 1 ) = (MT n (D)M~ 1 )Q' n . 

Hence MUM" 1 £ V{W') and T n {MDM~ 1 ) = MT n {D)M~ 1 . 

Corollary 2 . 2 . Let {i 7 'n}n>o be a pre-sequence of matrix orthogonal polynomials 
such that Fq is a C°°-function on ( a,b ). If F n D = A n F n for all n > 0, where D 
is a right handed linear ordinary matrix differential operator, then D = FqDFq 1 £ 
T>{FqWFq). In particular D is hypergeometric. Moreover if Q n = F n Ff x , then 
QnD = A nQn- 

Corollary 2 . 3 . Let {i 7 'n}n>o be a pre-sequence of matrix orthogonal polynomials 
such that Fq is a C°° -function on (a, b). If D is a second order differential operator 
whose leading coefficient is x(l — x) and F n D = A n F n for all n > 0, then 

(5) D = FqDFq 1 = d 2 x{\ -x)+ d{C - xU) - V, 

where C, U , V are constant matrices. 

If a pre-sequence of matrix orthogonal polynomials {F n } n >q comes from the 
irreducible spherical functions, of a fixed J\-type 7r, of a compact symmetric pair 
(G, K ) of rank one, then one can determine very explicitly a second order differential 
operator D = d 2 x(l — x) +dA\{x) + Aq(x) such that F n D = A n F n for all n > 0. In 
fact, this can be reached by considering the 7r-radial part of the Casimir operator of 
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G. For historical reasons let = Fq. Then 'FZ?'F 1 = D is equivalent to T being 
a solution of the following differential equations, 

(6) 2x(l — x)'&'(x) + ^{x)Ai{x) = (C — xU)^{x), 

(7) x{\ — x)^"{x) + T'(a;)Ai(j:) + T(a;)Ao(x) = -Ff(r). 

If one is able to compute T(a;) _1 , from © and (Q one can determine the constant 
coefficients C, U , V of the hypergeometric differential operator D given in (0 . 

If $ is an spherical function on G of A'-type 7r, then Q(kgk') = Tr(k)$(g)Tr(k'), 
for all g £ G and k,k' £ K. In particular if X £ t, then [Ad>](g) = k(X), 

for all g £ G. Let I(jr) denote the kernel of fr in D(K) and consider the left 
ideal D(G)I(n) and the quotient algebra D(n) = D(G) K /{D[G) K fl D(G)I(tt)). 
Then D(ir) is commutative, by hypothesis, because it can be embedded, using the 
Lepowsky antihomomorphism, into the commutative algebra D(A) <8> Eik1m(I4-)> 
M being the centralizer of A in K. 

Let D £ D(G). The 7r-radial part R(n,D) is a differential operator on A that 
satisfies R(tt, D)($\a) = D(d>)| J 4 . By taking an appropriate local coordinate x on 
A we can consider R(tt,D) as a differential operator on (0,1). If we conjugate 
R(ir,D) with T = Fo we obtain a differential operator D = ’FD’!/ -1 £ D(’I , I/F'F*), 
see Corollary 12.21 The map D H > D is an injective homomorphism of D(n) into 
Z^'FkF'l/*), which is not surjective because Z?(’FIF’F*) is not necessarily commuta¬ 
tive. An interesting challenge would be to find an extension of such homomorphism 
to a global algebra containing D(G) K in such a way to obtain an isomorphism onto 
D('FW'F*) (cf. H3j). 


3. An example 


The following example is taken from Theorem 3.2 in m- Let n £ No. Let us 
consider the following sequence of matrix polynomials 


F w (x) = 


3 F2 (_ iu ^+„ +3i2;x ) 


aFi ( ~ w ’ w 3 +n+3 \x) 
aFi'( “ 3 +n+4 ;*) 3A 2 (™+ c n + 3 ’ C+1 5 x) 


w is a nonnegative integer and c = ■ In particular 


Fo(x) = 


1 1 
1 1 — (n + 2)x 


We observe that det(F 0 (:r)) = — (n + 2)x. Hence F 0 (x) is singular only at x = 0. 

The rows of F w correspond to the irreducible spherical functions of type (n, 1) 
associated to the complex projective plane ^(C) = SU(3)/U(2) . 

From the orthogonality relations of the irreducible spherical functions, it follows 
that the sequence {-F,<,}tu>o is orthogonal with respect to the weight matrix 


W(x) 


(x(l — x) n+1 0 \ 

V 0 x(l-a :) n J’ 


supported in the closed interval [0,1]. 

The following theorem is obtained from Theorem 3.5 in [15J by taking a = n, 
j3 = 1 and using that F w (x ) = P w ( 1 — — x). 
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Theorem 3.1. The sequence {i 7 'iu(x)}u)>o satisfies the following three term recur¬ 
sion relation 


(1 x')F w (^x') — A w F w —.i(x') T B w F w [x ) T C U1 F W + i(x), 


with 


where 


A,,, — 


C w = 


By, = 


' w(w-\-n)(w-\-n-\- 2) 

(ii>+n+l)(2u>+n+2)(2'u;+n+3) 

o 


' (tti+l)(m+3)(tt;+n+3) 

(m+2) (2tt)+n+3) (2it;+n+4) 
_ m+3 _ 

, (u>+2)('u;+n+3)(2'it;+n+4) 

B w 

_ w-\-n -\-1_ 

(u>+l)(tt;+n+2)(2u?+n+4) 


_ w _\ 

(it;+l)(ii;+n+l)(2ii;+n+3) i 
w(w-\-2){w+n+l) I 

(it;+l)(2ii;+n+3)(2ii;+n+4) / 

0 'i 

( w+3)(w-hn+2)(w+n+4 ) I 
(m+n+3) (2u;+n+4) (2w+n+5) ) 

iu+n+3 _\ 

(-io+2)(ii;+n+2)(2it;+n+3) \ 

B 2 W 2 ) ’ 


oil _ (h>4~1) 2 (^+3) _ , _1_ 

(ir+2)(2ui+n+3) (2u>+n+4) ' (ii;+l)(ii;+2)(it;+n+l)(ii;+n+2) 

,_ (w+n)(w+n+ 2) 2 _ 

' (ty+n+l)(2i<;+n+2)(2ii;+n+3) ’ 

o22 _ _ (u>+l)(it;+3) 2 _,_ (m+n+l) 2 (m+n+3) 

(ii;+2)(2io+n+4)(2it;+n+5) ' (m+n+2)(2i<;+n+3)(2tt;+n+4) ’ 


Therefore the sequence {F„}iu>o is a pre-sequence of matrix orthogonal poly¬ 
nomials. Hence, according to Theorem 12.11 {Qy, = F u; F 0 " 1 } u ,>o is a sequence of 
matrix orthogonal polynomials with very nice properties which can be given explic¬ 
itly, since F 0 _1 can be computed. 

A careful computation leads to 


Q w (x) = F w (x)F 0 (x) 1 = 


Qf[(x) Ql> 2 {x) 
QSi(®) Ql’ 2 (x) 


where 


-w(w-\-n- 1-3) ( —'Uj+l, w-\-n +4 . ^ 


Qlli x ) ~ 3^2 V 3, l 
q ? 2 (*)= ~igy a) afi ( 

QSi («) = 2 Fi (-«'-™+"+ 4 
QM = ^J3F2 


—ih, iu+n+3, 2 . \ , w(i/i+n+3) 7-1 / — tu+1, tu+n+4 . 

>Xj + 3(n+2) 2^11 4 ’ X J, 

x) 


X - 


3 F 2 ( 


3 ’ “7 3(n+2) 

—w, w-\-n-\- 4, s«; + l . \ 

4, s« ) Xj , 


—w, w-\-n-\- 4, + 1 


;x), 


with s w = w(w + n + 4) + 3 (n + 2). It is worth to observe at this point that 
deg(F u ,) = w + 1 and that deg(Q u) ) = w. Moreover, the leading coefficients of F w 
and of Q w are, respectively, the following: 


0 0 \ 

LC(i^) = I g (-l) m+1 (m+n+3) m -|-i(c-|-ii)-|-l) I , 

V (3)m+ic / 


LC (Q w ) 


/ (~l) w (w+n+3) w 2 

I (2 -j-w)w! 

I ( — 1)™ (w-\-n-\-4) w w(w— 3) 
\ (3) m+1 (n+2) 


( — l) w (w-\-n+4) w (s w +w ) 
(n+2)(4)„ , 


In particular we see that the leading coefficients of F w are singular while those of 
Q w are nonsingular, for all w > 0. 
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The orthogonality matrix weight of {Q w } w >o is 


W' = F 0 WFq 


r(l -x) r 


2 — x 2 — (n + 3)x 

2 — (n + 3)x 1 — x + (1 — (n + 2)x) 2 


supported in the closed interval [0,1]. It is always interesting to wonder if a given 
matrix weight is irreducible or not. To this end we look at the commuting space C 
of W'. 

The commuting space C of W' is the set of all constant matrices such that 
TW'(x) = W'{x)T* for all x € [0,1]. If T = (*ij)i<»,j <2 belongs to C if and only if 
the following equations hold for all x, 

(1.1) : in(2 — x) + * 12(2 — (n + 3)x) = (2 - x)tu + (2 — (n + 3)x)*i 2 , 

(1, 2) : *n(2 — (n + 3)x) + * 12(1 — x + (1 — (n + 2)x) 2 ) 

= (2 - ar)* 2 i + (2 - (n + 3)x)t 22 , 

(2.1) : * 21(2 - x) + * 22(2 - (n + 3)x) 

= (2 — (n + 3)a;)*ii + (1 — x + (1 — (n + 2)a;) 2 )*i2, 

(2, 2) : *21 (2 — (n + 3)a:) + * 22(1 — x + (1 — (n + 2)x) 2 ) 

= (2 — (n + 3)x)*2i + (1 - x + (1 - (n + 2)x) 2 )i 2 2- 


From (1,2) it follows right away that *12 = 0. Then from (1,1) one obtains that 
*n = *n. Similarly, from (2,2) it follows right away that *22 = * 22 , and then that 
*21 = * 21 - At this point equation (1,2) has been reduced to 

*n (2 — (n + 3)x) = (2 — x)* 2 i + (2 — (n + 3)x)t 22 . 


Since this should hold for all x we obtain that *n = t 22 and that *21 = 0. Therefore 
C = R/, which is equivalent to W' being irreducible, see Corollary 1.13 of [22]. 


Since the irreducible spherical functions are eigenfunctions of the Casimir oper¬ 
ator of SU(3), it follows that (cf. Section 3 of [T5]) F W D = A W F W , w > 0, where 


D = d 2 x(l - x) + d 


2 — (n + 4)x 0 

0 2 — (n + 3)x 


+ If- 1 1 ~ a; 

1 1 —1 T x 


and 


A,„ — 


—w(w + n + 3) 0 

0 —ui(w + n + 4)^n — 2) ' 


Now we will compute FqDF 0 . 

FoDFo 1 

= d 2 x(l -x)+d(2 FqF 0 - 1 x( 1 -x) + F 0 


2 — (n + 4)a 


2 — (n + 3)s 


Fr 


-1 


F' 


(2 - (n + 4)x 

0 


0 




We have 




2 — (n + 3)x J 0 
— 1 fl — [n + 2)x —1 


-F 0 

x 


— 1 1 — x 

1 — 1 + X 


F~ l 
r o ■ 


(n + 2)a 


-1 


and FqF 0 = 


1/0 0 

-1 1 /’ 
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Fn 


2 - (n + A)x 0 \ ! 

2 — (n + 3)a; J 0 

1 (2n + 5 — (n + 2)(n + A)x 

1 — (n + 2)x 


n + 2 


-1 

2n + 3 — (n + 2)(n + 3)a; / 


F' 

n 


, /2 — (n + 4)x 


0 


0 


2 - (n + 3)xJ F ° x 


1 


0 0 
—2 + (n + 3)a; 2 — (n + 3)x 


Fn 


— 1 1 — a; 

1 — 1 + x 


F^ = 


0 0 
2 — (n + 3)x —2 + x 


This finally gives 


D = FqDFq 1 = d 2 x(l -x) + d{C - xU ) - V, 


where 


C = 


1 /2n + 5 -1 


+ 2 \—2n — 3 An+ 7 


, u = 


n + 4 0 

— 1 n + 5 


V = 


0 0 
0 n + 2 


Since F W D = A W F W and Q w = Fy^F^ 1 for all w > 0, it follows that Q W D = 
A w Qw f° r a ll w > 0. Let Qj tW (x) = (Qj.i(x), Qj t 2 ( 2 )) be the jtli-row of the matrix 
polynomial Q w , for j = 1,2. Then 


(8) x(l - x)(QU" + ( C* - xU^iQl J - (V* + \j, w )Qi w = 0, 

where Ai >u , = —w(w + n + 3) and \ 2 , w = —w(w + n + 4) — n — 2. 

The characteristic polynomial of C 4 is p( A) = (A—2)(A—4). Since the eigenvalues 
of C 4 are not in the set {0,-1,— 2 ,...}, the analytic solutions at x = 0 of (0) are 
given in terms of the hypergeometric function (cf. Definition 3 in |2I]) 

Ql w (x) = affr ( u *< yt c t Xi ' w !*) QiJ 0), 

with Qi iW (0) = 3 ( n ^ 2) (w+n+3,-w-n- 3), Q 2 , w { 0) = 3 ( n ^ 2 ) (~ , » ; -»-4, w+n+4). 
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